Effects of Chirped Laser Pulses on Nonclassical Correlation and Entanglement of 

Photon Pairs from Single Atom 
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We study the effects of arbitrary laser pulse excitations on quantum correlation, entanglement 
and the role of quantum noise. The transient quantities are computed exactly using a method 
that provides exact solutions of the Langevin field operators for photon pairs produced by a double 
Raman atom driven by laser pulses. Short pulses with appropriate chirping, delay and/or detuning 
can generate broadband photon pairs and yield results that provide insights on how to generate very 
large nonclassical correlation. We find that short pulses are not favorable for entanglement. The 
quantity was previously found to be phase-sensitive and this is used with the pulse area concept 
to explain the rapid variations of entanglement with pulse width and strength. Photon correlation 
and entanglement are favored by exclusively two different initial conditions. Analysis reinforces our 
understanding of the two nonclassical concepts. 

PACS numbers: 
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I. INTRODUCTION 

The intriguing nature of atoms and photons man- 
ifest through entanglement and quantum correlation. 
These nonclassical concepts carry subtly different phys- 
ical meanings although they are mathematically related 
for two- mode Gaussian state [1]. As pointed out by R. 
J. Glauber, quantum correlation exists when the den- 
sity matrix is not a product form Yl k |a^)(o^| (in co- 
herent state basis), rather in an incoherent superposi- 
tion J2k w k\ a k)( a k\ or mixed state [2]. Entanglement is 
commonly depicted as a superposition of two mutually 
inverted states. The physical mechanisms that create 
quantum correlation and entanglement may not be the 
same and not fully understood. 

We consider a single atom with four-level double A 
system (Fig. [1]) driven by laser pulses. This scheme is in- 
teresting for being able to generate controllable nonclas- 
sically correlated Stokes and ant i- Stokes photons from an 
atom [3] as well as extended medium Q . Previous works 
on single atom driven by continuous wave (c.w.) lasers 
were based on Schrodinger's equation 0], Q and master 
equation [6]. The former approach yields solutions only 
for c.w. lasers. The latter uses atomic correlation to ob- 
tain field correlations indirectly via quantum regression 
theorem [7[ . A general and more transparent method for 
computing correlations is based on Heisenberg-Langevin 
(HL) formulation, which can be extended to include spa- 
tial propagation. However, the previous studies have not 
considered excitations with laser pulses. 

Laser pulses provide extra degrees of freedom for co- 
herent control of nonclassicality. It is interesting to ex- 
plore the effects of pulse width (duration) and other pa- 
rameters such as chirp, phase and shape on quantum 
correlation and entanglement of photon pairs, as well as 
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laser pulses 



FIG. 1: (Color online) An atom with four- level double Raman 
scheme is driven by pump(p) and control(c) laser pulses, pro- 
ducing Stokes (k) and ant i- Stokes (q) photons with transient 
quantum correlation and entanglement. 



on the role of quantum noise [8]. Nonclassical photon 
pairs driven by laser pulses could be made sufficiently 
intense and serves as a new tool for doing quantum non- 
linear optics. A scheme to generate single-cycle photon 
pairs has been proposed [9[. For (short and intense) 
laser pulses, exact solutions cannot be obtained using 
the Schrodinger's equation. The quantum regression ap- 
proach may give numerical solutions for pulses, but can- 
not account for spatial propagation. The ultimate goal 
is to compute quantum correlation for extended medium 
driven by arbitrary laser pulses. 

As the first step towards this goal, we present a general 
method that gives exact transient solutions for quantum 
operators (a& and a q ) of photon pairs in HL equations 
for single atom. The solutions enable us to study, for 
the first time, the effects of various laser pulse parame- 
ters on the transient quantum correlation and entangle- 
ment of the photon pairs. We obtain interesting results 
that reflect the unique advantages of using pulses, thus 
showing the inherent advances in extending the previ- 
ous studies to pulsed excitations. We find that chirped 
and overlapping laser pulses can give an extremely large 
photon correlation. In addition, short pulses tend to in- 
crease the photon correlation but reduces the entangle- 
ment. We give physical explanations that shed insights 
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on the effects of finite spectral bandwidth on the non- 
classical properties. By analyzing the main results that 
contain rich underlying physics, we acquire conceptual 
insights on: a) coherent control of nonclassical photons, 
b) the differences between nonclassical correlation and 
entanglement of photon pairs and b) the effects of pulses 
on the role of quantum noise. 



II. THEORETICAL FORMULATION 



of operating parameters. We outline a method which pro- 
vides exact solutions for the field operators. 

First, note that the time evolution of the density 
matrix elements can be solved numerically as X(t) = 
(X(t)) = U(t)(X(0)), where X(0) = {cr af3 (0)} where 

aft = aa, ...a<i, 6a, bd, dd is the initial vector of the 

atomic operators and U(t) is the time-evolution matrix 
whose exact form can only be obtained numerically from 
the solutions p yx = (cr X y) of the density matrix equations. 
We then find the exact solutions 



The usual Hamiltonian for an atom in Fig. [T] interact- 
ing with radiation fields and two laser fields is 
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where fi p ( c ) is the Rabi frequency of the pump (control) 
field, Afc = Vk — ^abi A q = v q — co ac are the detunings 
of the Stokes (subscript-k) and anti-Stokes (subscript- 
q) frequencies, g* are the usual atom- field coupling 
coefficients. 

It is straightforward to obtain a set of sixteen HL equa- 
tions including noise operators as represented symboli- 
cally by 



d_ 

dt 



—+°xy(t) 



r xy&xy (t) 
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+ F xy {t) (2) 



where F T 



Fxy& 



iA xy stands 



xy J - xy^ ' xy and r ^ xy Ifxy 

for the complex depopulation or clecoherence rates, C, 
is the coefficient that depends on the laser fields and F xy 
represents the noise operators. The HL equations rep- 
resented by Eq. ([2j) are obtained after eliminating the 
solutions of the field operators dk and d q in the Heisen- 
berg equation -^O = ^ [O, H] for atomic operators using 



d 
dt 



(3) 



4w = 4w 



K m (t)X m (0) + / K m (t')F m (t - t')dt' 
Jo 



(5) 
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Qn(t)X n (0)+ Q n (t')F n (t-t')dt' 
Jo 



(6) 



where 
rt 



K m (t) = fiU^mWdt and Q n (t) = 

Jq U9 y7l (t f )dt f . This technique for obtaining exact opera- 
tor solutions is applicable to an arbitrary time dependent 
laser fields interacting with an atom. It may be extended 
to molecules in future works. 

The single operator expectation values in thermal 
vacuum are = -ig% J2m=i K m (t)(X m (0)) and 

(d q (t)) = ig q Y^m=iQm(t){X m (0)). From the solutions 
Eqs. ([5]) and (|6|), we can compute the expectation values 
for the products of two operators. For example, 
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(a q (t)a k (t)) = -g k9q £ [Q m (t)K(t)(X m (0)Xl(0)) 

m,n=l 
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Q m (T - t')K n (t - t')2D™ n (t')dt'} (7) 



d_ 

dt 



°q(t) = ig q Vca{t) 



(4) 



where crfe(t) = due %Vkt and d q (t) = d q e 



III. METHOD FOR EXACT SOLUTIONS 
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n k =nt + \g k \ 2 £ [K m (t)K* n (t)(X m (0)Xl(0)) 
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K m (t-t')K* n (t-t')2D™ n (t')dt'} 



(8) 



Exact analytical solutions for dk and d q can be found 
from the exact solutions of (J a b(f) and & ca (t) (since these 
equations are linear) but the solutions are very cumber- 
some and not helpful. Adiabatic approximation may sim- 
plify the solutions but it will considerably limit the range 



where = (d\(t)dk(t)) is the Stokes photon num- 
ber, (4(0)a fe (0)) = nf and (aj(0)a g (0)) = n\ h are the 
mean thermal photon numbers at initial time. The nor- 
mal ordered and antinormal ordered diffusion coefficients 
are defined as £>»„(*) = (F%(t)F n (t)) and D»» (t) = 
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(F m (i)F^(i)) respectively. The final terms that con- 
tain the diffusion coefficients are the noise (superscript 
V) parts, i.e. n k and n™, while the former terms are 
referred as the boundary parts n b k and n h q . We also 

use (JafiCJ^ = ((J-ypd'pay = ^a-f- We have assumed 
that initially the modes are statistically independent: 
(a q (0)a k (0)) = (a q (0)al,(0)) = and are not squeezed, 
i.e. <a2(0)> = (a g (0)> = 0. 

IV. NONCLASSICAL QUANTITIES 

The average value of the paired operators are com- 
puted to obtain mean photon numbers, photon-photon 
correlation and to determine entanglement. The effects 
of finite pulse duration, chirping and pulse sequence on 
these quantities can be studied. We also analyze the ef- 
fects of pulses on the contributions of the noise operators 
compared to the initial operators. 




a) c.w. lasers b) pulsed lasers 



FIG. 2: (Color online) Effect of pulses on the Cauchy-Schwarz 
correlation for DRR (A p — A c = 0,Q P — Q c ) scheme with: a) 
c. w. lasers, b) identical and coincident laser pulses of width 
a c = l/15r ac (profile shown by the dashed line). Note the 
substantial increase in the nonclassical correlation (g cs > 1) 
for the case driven by short pulses. Parameters used are 
Q p , c = 10r ac , no decoherence 7& c = and the initial pop- 
ulations /?cc(0) = pbb(0) = 0.5. We assume Gaussian profile 
exp[— (t — to) 2 /a^] for the pulses. 



A. Nonclassical Photon Correlation 

The solutions Eqs. ([5j) and ([6]) are linear combina- 
tions of initial field operators and atomic noise operators. 
Thus, the two-photon correlation function can be decor- 
related into terms with paired operators. The presence 
of nonclassical two-photon correlation between the Stokes 
and ant i- Stokes photons is determined by the condition 
gCS y 2 for the Cauchy-Schwarz correlation 

(a{(t)al(t)a q (t)a k (t)) 
g cs (t) = ± -L 

I (a g afc) | 2 + | (o$Qfc) | 2 + (nfc) (h q ) 
^{Kal}P + 2(n fe ) 2 }{|(d2)|2 + 2(n g ) 2 } 

where G^p (t) = (^d^(t)d^(t)dj(t)dj(t)^ (j — k,q) and we 

omit the t dependence to simplify the notations. Note 
that Eq. (|9]) corresponds to transient coincident (joint) 
two-photon detection (with no delay between the Stokes 
and anti- Stokes) which is mathematically related to an 
entanglement criteria Eq. JTDJ) below. 

B. Entanglement 

The above solutions Eqs. ([5]) and (|6]) are also used 
to identify the presence of entanglement via D(t) = 
((Au) 2 ) + ((Av) 2 ) < 2 [11]. This is a sufficient and 
necessary to criteria for two- mode Gaussian states. We 
compute the paired operators: (d q {t)dk(t)), (d\,(t)dk(t)), 
{a\(t)a q {t)), (a q {t)a\(t)), (d 2 k (t)), (a 2 q (t)) just like Eqs. 
© and J5J). Note that the finite mean thermal photon 
numbers nj^ e or n q he effectively increases the value D. 




a) Variation with pulse width b) Variation with control field 



FIG. 3: (Color online) Variation of the correlation profile 
with: a) pulse width a c (= <J P ) for resonant case A p , c = 0, 
with Q P;C = 10r ac and b) Rabi frequency Q c (= Q p ) for de- 
tuning A P;C = — 100r ac and width a c = l/15T ac - Other 
parameters are the same as in Fig. [2] 



Thus, entanglement may not be obtained when there is 
a significant number of thermal photons. 

Let use recall the steady state for two-photon state 
1$) = Ek,q C kq|c,l k ,lq) H ■ For Raman-EIT scheme 
which shows nonclassical correlation, the coefficient is 
inseparable Ck q ^ C^C^ [5]. Thus, one might expect 
to find entanglement in the steady state as well as in 
the transient. However, we find no entanglement for the 
Raman-EIT scheme at all times and with any control field 
Q c and phase (j) c of the laser. We emphasize that the in- 
separability of Ckq only implies correlation, and may not 
give entanglement depending on the complex phase <\> of 
(a q (t)a k (t)) through Q 



D = 2[l + n k +n q + 2^n k n q (gW - l) C os <j> kq (10) 

which relates the quantity D with normalized correlation 

# (2) = (a\a\a q a k )/n q n k . 
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V. RESULTS 

Based on the results, we discuss the effects of laser 
pulses on nonclassical photon correlation, entanglement 
and the contribution of quantum noise. 



A. Nonclassical photon correlation 

Several interesting effects are found for excitations 
with ultrashort pulses. First, let us compare the results 
(Fig. [2]) with c.w. lasers and pulsed lasers for identical 
and resonant lasers: A p = A c = 0,Q P = Q c referred 
to as double resonant Raman (DRR). Continuous wave 
lasers give a Fresnel-like oscillations between classical and 
quantum regimes, showing no apparent nonclassical cor- 
relation (Fig. Ek). The oscillations have a period 7r/^ P)C , 
are due to optical nutation [TO] which do not show up in 
the case of pulsed excitations due to the adiabaticity. 



1. Effects of finite pulse width 

Laser pulses give nonclassical correlation much larger 
(Fig. [2)3) than the correlation from c.w. laser. This is 
one of the main results. The correlation increases sub- 
stantially for pulse width a c shorter than 0.1 /T ac (Fig. 
[3^i) . The physical explanation comes from two points: 
a) A short pulse carries a large bandwidth, b) Correla- 
tion is a measure of the amount of spectral interrelation 
between the Stokes and anti- Stokes photons. 



2. Effects of pulse detuning 

Pulses with finite detunings A PyC » Q PyC produce 
large nonclassical correlation (Fig. [3)3). Here, the band- 
width of the generated pulses of photons are enhanced by 
another mechanism (in addition to the short pulse dura- 
tion); i.e. the spontaneous Raman processes from the 
two laser pulses. However, for c.w. case, the correlation 
does not depend on the Rabi frequencies regardless of 
whether there is detuning or not. 



Effects of Chirping 




a) Spectral effect due to chirping b) a p =a c , coincident pulses 




c) a p = -a c , coincident pulses d) a p = -a c , overlapping pulses 



FIG. 4: (Color online) a) Spectral content of chirped carrier 
frequency, computed from the Fourier transform of exp[i(zy c £+ 
a c t 2 )]. Effects of chirping on correlation profile for: b) co- 
incident pulses and identical chirp, c) coincident pulses and 
opposite chirp, d) partially overlapping pulses and opposite 
chirp. Other parameters are the same as in Fig. [2J). 

photons that resulted in a large correlation. The corre- 
lation can be further enhanced by displacing the con- 
trol pulse slightly such that both pulses overlap par- 
tially. This yield an extremely large transient "corre- 
lation pulse". If the pulses are entirely separated, the 
chirping has no enhancement effect on the correlation. 

The Raman-EIT scheme has a very large g cs at times 
much smaller than the lifetime but decays rapidly to 
zero for large times, corresponding to steady state photon 
antibunching found in previous studies Q. 



B. Entanglement 

Figure [5] shows the entanglement for different pulse 
sequences for DRR scheme. 



3. Effects of chirping and delay of pulses 

We now analyze the effects of chirping of the pulses. 
We consider the chirped carrier frequency of the form 
v x (t) = v x + a x t (x = p,c). Figured^ serves to remind 
the spectral content of positively chirped and negatively 
chirped pulses. When the pulses are chirped identically 
(a p = a c ), there is no significant effect on the correlation 
(Fig. SJd). For opposite chirping (a p — —a c ) and suffi- 
ciently large |ctp, c |, the correlation shows a peak. This 
can be understood qualitatively from Fig. [4^. The two- 
spectral "wings" span an effectively broad spectrum of 



1. Favorable initial condition 

The symmetrical initial condition p cc (0) = pbb(0) — 
0.5 that gives quantum correlation does not give entan- 
glement. We obtain entanglement with a different initial 
condition p cc (0) = 1. This initial condition ensures that 
Stokes photon is produced first before an anti- Stokes pho- 
ton since there is little or no population in level b. But 
a more important point that explains entanglement is 
that, the state of the atom which generates the anti- 
Stokes photon carries quantum information about the 
Stokes photon. The simultaneous excitations by lasers 
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counter-intuitive sequence 

intuitive sequence 




FIG. 5: (Color online) Entanglement between the Stokes and 
anti- Stokes photon pairs for the DRR scheme using: a) con- 
tinuous wave (c.w.) lasers, b) coincident pulses, overlap- 
ping pulses with c) intuitive sequence, d) counter-intuitive se- 
quence. The overlapping pulses are shown on top of the figure. 



Parameters are: cr p , c = 
7 &c — and p cc (0) — 1. 



Q = 3or. 



(2.ir oc 



o,a, 



3or Q 





a) Variation with pulse width b) Variation with control field 

FIG. 6: (Color online) Variation of the entanglement pa- 
rameter D for DRR scheme with: a) pulse width cr p , c with 
30r ac and b) Rabi frequency Q pc with pulse width 
(2.ir ac ) _1 for coincident pulses centered at 
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c) spontaneous Raman 



1 2 tr 3 

ac 

d) spontaneous Raman-EIT 



FIG. 7: (Color online) Mean number of Stokes photon com- 
puted from the solutions with boundary (initial) operators 
only (n|), noise operators only (njj) and both boundary and 
noise operators (rife) for: a) resonant (A p = 0) Raman with 
c.w. pump (O p = 5r ac ), b) resonant Raman with pulsed 
pump, c) spontaneous Raman (A p = — 50T ac , Q p = 5r ac ), 
b) spontaneous Raman-EIT (add a resonant control laser, 
C = 10r ac ). We assume jbc = and the initial population 

Pcc(O) = 1. 



0.5/r ac - Other parameters are the same as in Fig. [2] 



Figure [6] illustrates how the entanglement parameter D 
varies with pulse width a PyC and laser field Q PyC for DRR 
scheme. The oscillatory feature implies that the presence 
of entanglement could be very sensitive to the variations 
of <j PiC and fl PjC - Our analysis shows that the oscillations 
in the transient entanglement can be understood via the 
ulse area concept. The parameter D is phase sensitive 
l| and would be affected by the amplitude cr PiC ft PiC of 
the pulse area which manifests as a phase. 



C. Quantum noise 



ensures that both Stokes and anti- Stokes are driven with 
the same phase. This explains how identical lasers con- 
figuration and Pcc(O) = 1 can generate entangled photon 
pairs. 



2. Effect of pulses on entanglement 

The increased spectral bandwidth due to the pulses 
may be favorable for the correlation (see Fig. [2)3) but 
not for the entanglement. Coincident pulses give entan- 
glement while partially overlapping pulses in both intu- 
itive and counter-intuitive sequences do not give entan- 
glement. 



Figures [7ti,c show the Stokes photon number for single 
Raman transition with only the pump field. As expected, 
the contribution of quantum noise is significant for reso- 
nant pump field (Fig. [7^) but negligible for spontaneous 
Raman (weak and far detuned pump) (Fig. Et). This 
shows that quantum noise is important and contributes 
to the photon numbers when the field is resonant and 
gives rise to spontaneous emission. 

If a pump pulse is used instead, the Stokes photon 
number increases significantly but the proportion of the 
noise contribution is reduced (Fig. 03). This can be 
understood: the pump pulse corresponds to a range of 
frequencies, but only the resonant frequency contributes 
significantly to the noise through spontaneous emission 
while most frequencies that are off-resonant undergo 
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spontaneous Raman process which gives very little noise. 

When a resonant c.w. control field is applied (in ad- 
dition to the pump field) we have ant i- Stokes photons in 
the Raman-EIT (electromagnetic induced transparency) 
scheme (A p » f2 p , A c = 0) [4]. The contribution of 
the noise to the Stokes photon number becomes signifi- 
cant despite the large pump detuning (see Fig. [7]i). This 
seems surprising, but can be explained as due to the noise 
in population level c due to the spontaneous emission 
noise of the anti-Stokes photons. 

VI. CONCLUSIONS 

To summarize, we have presented a method to obtain 
exact solutions of photon operators that include quan- 
tum noise. The method enables exact computation of 
a variety of quantities that involve products of photon 
operators such as photon correlation, entanglement and 
squeezing in the transient regime. The results enable us 
to study the effects of pulsed excitations on the nonclas- 
sical photon correlation and entanglement. Short laser 
pulses produce much larger quantum correlation com- 



pared to c.w. lasers. The correlation profile is essentially 
independent of laser field for resonant pulses. Finite de- 
tuning of the pulses can very large correlation. Most 
importantly, we have studied how the laser pulses pa- 
rameters like chirping, duration and pulse sequence can 
give an extremely large nonclassical photon correlation. 
Chirping of the pump and the control lasers in an oppo- 
site manner increases the correlation. The correlation 
increases further for partially overlapping pulses. We 
provided explanation based on the spectral content and 
large bandwidth of the laser pulses that conspire to en- 
hance the photon correlation. We learned that the initial 
condition of the atom determines the transient entangle- 
ment and correlation. We explained how identical lasers 
configuration and initial condition p cc (0) = 1 can gen- 
erate entangled photon pairs. We apply the pulse area 
concept to explain why the transient entanglement varies 
(oscillates) with pulse width and the laser field. We have 
also found that the laser pulses tend to reduce the con- 
tribution of quantum noise to the photon number. All 
these results provide conceptual insights that show the 
significance of pulsed excitations on quantum properties. 
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